(N 

C/3 



S 



Testing a Variational Approach to Random 
Directed Polymers 

Andrea Pagnani* 

Dipartimento di Fisica, Universita di Roma Tor Vergata 
Via della Ricerca Scientifica 1, 00133 Roma (Italy) 



t^ : 1st February 2008 

fj ■ Abstract 

The one dimensional direct polymer in random media model is 
investigated using a variational approach in the replica space. We 
QQ . demonstrate numerically that the stable point is a maximum and the 

cn I corresponding statistical properties are for the delta correlated poten- 

^T^ ' tial in good agreements with the known analytic solution. In the case 

f~^ ■ of power-law correlated potential two regimes are recovered: a Flory 

0^ i scaling dependent on the exponent of the correlations, and a short 

- ' range regime in analogy with the delta-correlated potential case. 



g. 

I ■ 1 Introduction 

Q , A challenging open problem in statistical mechanics is the behavior of fluc- 

f-^ I tuating random manifolds with quenched disorder (for a review see JIl). New 

insight has been obtained in the framework of replica formalism and, at 
^ ■ least at the level of the tree approximation (Hartree-Fock) which is the ex- 

act solution in the limit of infinitely many dimensions, the problem has been 
hopefully solved [|^ . Nevertheless the approximation turned out to be maybe 
too drastic for take in account the behavior in all dimension and for different 
structures of correlation of the quenched noise. For instance the prediction 
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for the wandering exponent ( of one- dimensional direct polymers in random 
media (DPRM) are wrong even in the case of short range correlations of the 
random potential. 

In the following we will restrict to DPRM, which due to their one-dimensional 
structure are simpler then random manifolds with generic dimension. 

After this short introduction we will present some details on the model 
and we will define some basic mathematical tools in order to introduce in 
section 3 the variational approach developed in and reformulating it in 
a suitable form to implement a numerical simulation. Section 4 is divided 
in two parts: in the first one we check the agreement of our results with 
Kardar's analytic solution [Q] for the delta-correlated potential and in the 
second part we study the power-law correlated potential for different values of 
its exponent. Two different scaling regimes are recovered: a short-range one 
analogous to the delta-correlated potential, and a Flory like regime in wich 
the wandering exponent depends from the power-law correlation exponent in 
agreement with the one loop renomalization group analysis presented in [Q. 
Finally in section 5 we present comments and perspectives. 

2 The Model 

Knowledge of end point probability distribution in DPRM is of fundamental 
importance. In principle it is always possible to determine this probability 
distribution Pt[^] starting from the explicit structure of the Hamiltonian of 
the system which depends on the space coordinate set x, the time t, and the 
random potential rj: 



«M = /''' {if + "('=•*') « 

The related partition function satisfies the imaginary time Schrodinger equa- 
tion: 

^^|l^ = lv^Z(x,t) + ry(x,t)Z(x,t) (2) 

The polymer's end point probability distribution function p(x, t) is readily 
obtained by: 



We want to implement a variational approach in order to determine a func- 
tional of the probability Pt [p(x, t)] which allows direct calculation of the 
quenched noise average of a generic observable O by the functional integra- 
tion: 

0[p(x,t)] = J rf[p]P4p(x,t)]0[p(x,t)] (4) 

Quenched averaging is believed to be a rather difficult problem in spin-glass 
and DPRM models. Replica technique consists on calculating quenched 
average of the replicated partition function. We can therefore introduce 
Z({xa}, t) = 2'(xi, t)Z(x2, t)...Z(xn, t) as the solution of the imaginary time 
Schrodinger equation: 



dt 



^ j=l a=l ^•^i,a a<b 



Zi{^a},t) (5) 



where V^(x — y) = ?7(x, t)r7(y, f). A couple of explicit solutions to (|) are 
already available for V{x — y) = 5{x — y) P and V{x — y) = {x — yY 0, 
both of them in the one dimensional case d = 1. Unfortunately those are 
the only known analytical solutions of the problem which in general can be 
only studied by mean of renormalization group approximate technique 0. 
A particular solution of (||) can be decomposed over the eigenstate's set of 
the Hamiltonian: 

Z({x4,t) = ^QV[/,({x°4)v[/,({x4)e-^'* (6) 

i 

From the previous expression it is evident that all the excited states decrease 
sub-exponentially compared with the ground state, and then the thermo- 
dynamic behavior of the model will be ruled in the asymptotic regime by 
the structure of the ground state. In this regime is possible to calculate for 
instance: 

p(xi)p(x2) = / dP[p] p(xi)p(x2) OC / dX3dX4...(iXn V^(xi,X2,...,Xn) (7) 

where '?/'(xi,X2, ...,Xn) is the ground state of (||) and dP[p] = lirat^oo dPt[p]. 

Since we know the exact solution of the one dimensional delta-correlated 

potential model, a careful analysis of its ground state structure should be 

illuminating for further development of the variational principle. In fact a 



straight application of the Bethe ansatz - which turns out to be the exact 
solution of the delta-correlated model - shows that the ground eigenstate is: 

-A\x. 



with energy: 



i){Xi,X2,..-,Xn) 



^W 



a<h 



Xh\ 



n[n 



l)A^ 



(8) 



(9) 



where A is a constant depending on the parameters of the Hamiltonian. So, 
following Mezard et al. in [Q , in the general case a reasonable choice should 
be restrict the space of n-particles test functions to the following: 



ll){Xi, X2, ...,Xn) = Y[ fi^a - Xb) 
a<b 



(10) 



A direct calculation shows that this is tantamount to choosing the following 
Gaussian probability distribution Pg[0] for 0(x) = ln(p(x)): 

1 



Pr. 



oc exp 



rf^a;rf^0(x)^(x-y)0(y) 



(11) 



In fact performing a Gaussian integration we obtain: 



(i[0] exp 



La=l 



^|rf'^xA0(x)(7(x-y)0(y) 



a,b 



Yb) 



(12) 
and the following relation holds: 

I d'z ^(x - z) ln[/(z - y)] = 5\^ - y) (13) 

In the previous one dimensional case g must be chosen in such a way that: 



dxdy (t){x)g{x - y)(t>{y) 



dx-\V<P\' 



(14) 



It's now reasonable to guess for our ground state a functional structure of 
the following type: 



^(X1,X2, ...,X„) OC / dR[fj] /i(xi)/i(x2)...yU(Xn) 

The distribution functions P[p] and R[p] are related by: 

/i(Xi)/i(x2) 



p(xi)p(x2) = / dP[p] p(xi)p(x2) 
where /[p] = / d'^x p(x 



dR[ij] 



m' 



(15) 



(16) 



3 The Variational Principle 

Since now we have recast our problem in a (hopefully) suitable form to im- 
plement the variational principle which consists on minimizing the quadratic 
form {il)\H\%l)) , where H is the replicated Hamiltonian operator which appears 
in the square brackets of @. In order to obtain normalized eigenstates we 
have to enforce this condition by the introduction of a Lagrangian multiplier. 
Following this procedure described in |^ we easily obtain: 

{il:\H\%l)) = f dR[iJ.]dR[a]nn[l^,(T] (17) 



where 



H.[/.,a] = ^/rf'^xEg;^/[H"-^+ (18) 

n(n -1)1 d''xd% ^(x)^(y)l^(x - y)a(x)a(y) J[H""' 
The normalization in terms of our trial functions is: 

(^1^) = f dR[ij,]dR[a] J[H" (19) 

In this new representation the Schrodinger equation reads: 

/ dR[a] Hnil^, a]=En f dR[a] /[H" (20) 

It's now possible to consider the n —>■ limit of (^). In this limit we 
can eventually implement the variational principle in terms of the stationary 
point of the free energy functional F[R] = J dR[fi]dR[(T]T-C[fi,a], and TC is 
now defined by: 

+ j d-^xd^y /i(x)/i(y)a(x)a(y)V^(x - y)/[H"' (21) 

As far as replica symmetric approximation is concerned, the problem is recast 
just in finding the proper (^-function - as defined in (pJ]) - which is stationary 
against variations of the following functional: 

y'rfP,[0i]rfP,[</)2]H[exp(</)i),exp(02)] (22) 



Performing analytic calculation on this functional is hopeless, so we are forced 
to try a numerical approach in order to calculate the critical exponents of 
the model. 



4 Numerical Analysis 

4.1 The delta-correlated potential 

As already mentioned, 1-dimensional DPRM with delta-correlated random 
potential is explicitly solvable and so it is a good testing ground for our 
variational approach. For technical reasons it is simpler to consider the case 
in which a lattice spacing is introduced. The Hamiltonian, taking lattice 
spacing, a = 1 is: 

1 ^ 
^ i=l 

+ G^fiifijaiajVijI[fi(j]~'^ (23) 

id 

where in the delta-correlated case Vij = 5ij. The continuum limit can be 
reached sending the lattice spacing to zero or sending the coupling constant 
measured in lattice units to zero. Indeed the dimensionless coupling constant 
is Ga^. Although the wandering exponent is supposed to be independent from 
the lattice spacing, the continuum limit may be interesting as a check to our 
computation. 

In order to implement a numerical integration of (^2]) we need as a first 
step to generate the fields which follow the distribution Pg [(f)] . The Fourier 
representation of this probability distribution is: 

P,[0]rf[0] ex n d[(l)k] exp{gk\(j)k\^) (24) 

fc<A 

where k G {{iTn) / L}n=i^2,...N , N are equally spaced intervals in which we 
consider discretized the field </'(x) x G {0, L} and (pk are Fourier coefficients 
of 0. The integration in Fourier space "decouples", in the sense that now 
the probability distribution is the product of the individual /c-wave vector 
probability distributions. 



Assuming that the wave function on the lattice is the same as in the 
continuum (the difference among the two going to zero in the continuum 
hmit), a httle algebra reveals that: 

Y[,-AK-x,\ =Jd[<p]f[expU{x^)-Jdx ^^1 (25) 

a<b a=l (. ) 



which means that for delta-correlated noise we straightforwardly obtain: 

-' 11 r — 

A:<A 



dP[<P] oc exp [jdx ^^1 d[<P] oc n exp {^^| ^W (26) 



From the structure of (|26|) we can reasonable guess the following structure 
for the probability distribution: 

dP[0] oc n exp j-^^!^| d[<t^,] (27) 

fc<A I ^ J 

in which both a and A are to be considered as variational parameters related 
with the fluctuations of d), since: 



((0(x) - <P{y)Y) = J dPg[<pmx) - <P{y)Y = \\x - yr (28) 

Comparing (p6| ) with ( pTf ) it's easy to deduce that the delta-correlated case 
is obtained setting a = 1, so we used this result to a first tuning of the 
variational parameters. We have generated statistically uncorrelated random 



numbers distributed as (27) with a = 1. Back Fourier transforming via 
FFT algorithms [Q the coefficients (pk, we are able to obtain statistically 
uncorrelated realizations of 4>{x). Let F be the set of M realizations of 4>{x), 
we have found a suitable interval of values for A in which: 

^ E (0((^a) - <P{ya)f = C\x-y\ C = 0(1) (29) 

asr 

The 0(a;) fields built following this prescription are the configurations used 
for calculating the value of the free energy functional Ti defined in (^). It is 
obvious that such a value should depend on the realization of the fields and 
so we have to average over a significative number of this realizations, or more 
formally, |0a) and j-^a) being the a-th of M realizations, we can compute: 

1 M 

{P) = TpT. i^aini^b) (30) 

a, 6=1 
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Figure 1: Fluctuations of the field vs. distance averaged over 10^ realiza- 
tions of disorder for (a = 1.0, A = 1.0). 



where in this case we use the delta-correlated potential V {x — y) = 6{x — y) . 
It's easier to use intensive variational parameters, so we redefine (^7D as: 



dPM 



n 

fc<A 



fcjexp 



-iV3/2 



Xk 



2a' 



(31) 



where a and A are new variational parameters. 

A crucial question now is whether the free energy functional should be 
maximized or minimized, since from a strictly formal point of view we are 
considering the n —>■ limit, a regime in which thermodynamic behavior is 
not evident at all. In the general theory proposed in it is proved that 
Hartree Fock approximation consists in a saddle point on the number of 
components of the field going to infinity, and the solution is a maximum 
of the functional. A possible argument which explains this rather unusual 
feature of replica theory is that the Hamiltonian (j^) is translational invariant, 
so if we choose the center of mass reference frame we are left with n — 1 degree 
of freedom (a number which is negative as soon as n < 1 !). This is the regime 
in which the stationary point becomes a maximum of the functional. 

The numerical simulations have clearly pointed out that the stationary 
point of the functional is really a maximum confirming the reasonment. Let 
us first present what happens for the coupling constant G = 1. Setting 
A = 1 the free energy functional develops a maximum for a = 1.1. The 
maximum value is obtained by varying also the A parameter. Iterating this 
procedure in the interval A G {0.8 : 1.3} we find that as the value of A grows 
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Figure 2: Plot of free energy F vs. a. The maximum value of F is obtained 
for (a = 1.0, A = 1.4). 



inside the interval, the free energy maximum shifts from (a = 1.15, A = 1.0) 
to (a = 1.0, A = 1.4) increasing monotonically its value. Over the point 
(a = 1.0, A = 1.4) we observe a decrease of functional. The range of the 
free energy functional in this window of parameters belongs to the interval 
{0.24 : 0.26}. 

The value given in (P) is F = 1/3 which is of the same magnitude of our 
simulation but significatively outside the error bar. This difference from the 
analytic value must be addressed to the fact that we are replacing a contin- 
uous model with a unitary step discretization so that only wave functions 
which vary on a length scale bigger than our discretization make sense. It 
is therefore reasonable to introduce a modulation on the coupling between 
replicas in (^I]) tuning the coupling constant G. The corresponding ground 
state energy is now G^/3. As far as values of G < 1 are concerned, the 
binding between replicas is lowered, and bound states with variations on a 
wider length scale are favored. We can't consider anyway too small values of 
G since in this case it's impossible to build a bound state and the stationary 
wave function turns out trivially to be constant. 

Following the previous considerations for each value of the coupling con- 
stant G in the interval {0.5 : 1.2} we have generated 10^ realizations of cf) with 
(A = 1.4, a = 1.0). As showed in figure (Q) we were able extrapolate from 
the linear fit of E{G)/G'^ vs. G^ the value of the abscissa F = 0.331 T 0.005 
which fits surprisingly well the analytic value of 1/3. 

The variance around mean value overestimates the statistical error related 




Figure 3: Plot of E{G)/G^ vs. the square of the couphng costant G^. 
Here A^ = 64 and (A = 1.4, a = 1.0). From the hnear fit extrapolate 
limG2^o E{G)/G^ = 0.331 T 0.005. 



with our numerical estimates of free energy, since as we can see from the 
normalized histogram of energy distribution in figure, relevant power-law 
tails are present. In order to give a reasonable estimate of statistical errors, 
we have generated for each experimental point 10^ realizations of free energy, 
which we have randomly divided into 10 subsets of 10^ elements. If - say K 
- is the number of subsets, the statistical estimate of variance made using 
the whole number of data we have, is i^ — 1 times bigger than the estimate 
given from averaging every variance in the subsets built under the previous 
prescription. 

It's worth noting that the variational parameter a is strictly related with 
critical exponents of the model like for instance the wandering exponent C 
defined by the relation (|x|) ~ t''. As pointed out in [§ 0(x) ~ x'^/t, and since 
from ( pSD we know that 0(x)^ ~ x", is possible to argue that ( = 2/(4 — a). 
In one dimension delta-correlated potential we found the stationary point for 
a = 1, which means ( = 2/3, in complete agreement with Kardar's analytical 
solution 0. 

4.2 The power-law correlated potential 

The verification test implemented for delta-correlated potential of the vari- 
ational principle is in good agreement with theoretical predictions. Unfor- 
tunately no solution is available for more general correlation of the random 
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Figure 4: Normalized histogram of the free energy frequency distribution for 
different G obtained by 10^ iterations of (H). Here N = 64, (a = 1.0, A = 
1.4). 



potential like power-law. It would be interesting for instance to know how 
the wandering exponent depends on power-law correlation's exponent. An 
approximate answer to this question arise from the one loop renormalization 
group analysis presented in |^ , which lead to the following scenario (defining 

{V{x)V{y)) (x\x- y|-(2p-i) = |x - y\-''): 

- V > 1/2 Renormalization group flux with one stable point which 
scales to the delta-correlated potential behavior. This regime is char- 
acterized by a fundamental length scale bigger then the length scale 
over which the effects of the correlation of the potential are sensible. 
Here C, = 2/3 {a = 1) independently from u. 

- z/ < 1/2 In this regime correlations extends on the same length 
scale of correlations and wandering exponent depends on p. It is known 
sometimes as the Flory's regime and the relation found is C(^) = 3/(4 + 
z/) which following our notation also reads a(z^) = (4 — 2z/)/3. 

In order to verify the approximated predictions above, we have inserted in 
the functional (^) a potentential of the following type: 



V{x-y) 



\x-y\ 



V 



V 



2p 



(32) 



It has to be notice that functional structure of the correlations has been 
chosen in order to obtain a logarithmic behavior in the limit z/ — ;> 0, so that 
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Figure 5: On the left is prented the behavior of parameter a which maximizes 
the functional vs. u {N = 64). A transition to the Flory long-range scaling 
regime is clearly observable for u > 1/2. On the left a zoom on the transition 
point is presented for A^ = 128, 256, 512. We can observe a more sharp 
transition from short to long range behavior as the volume increases. 



formation of bound states is still allowed in this limit. 

The interval u G {—2 : 10} has been subjected to investigation. Since 
negative values of the exponent u obviously correspond to correlations which 
grow with the distance, the finite size effects become more and more relevant. 
That's why was impossible to analyze the behavior of the system for values 
of z/ < -2. 

For each value of z/ we have calculated the corresponding pair {a, A) which 
maximize (^), calculating free energy and errors with the probabilistic pre- 
scription stated above. The results are presented in figure (for details see 
caption). From figure (|^) we find a rather smooth behavior around the tran- 
sition point u = 1/2. Iterating our analysis for polymers of different lengths 
(A^ = 128, 256, 512) we can observe a more sharp transition since finite size 
effects are less important. Figure (^ shows that for larger volume the dif- 
ferent curves develop a real discontinuity in the first derivative in i/ = 1/2 in 
agreement with the approximate results of Q. 
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5 Conclusions 

The variational method developed seems to correctly describe unidimensional 
polymers in random environment. For delta-correlated potential, stationary 
point is effectively a maximum of our replicated free energy functional in the 
?T, — > limit according with replica theory. The value of the maximum free 
energy is in very good agreement with Kardar's solution [^. 

The power-law correlated potential has been subjected to investigations 
and the results are in agreement with the perturbative result in |^: for 
z/ > 1/2 the wandering exponent equals the delta-correlated one (C = 2/3), 
and for v < 1/2 the standard Flory regime is recovered {({i') = 3/(4 + u)). 

A number of recently published papers pointed out that one dimensional 
DPRM are characterized by replica symmetry breaking. Our formulation 
can be in principle enriched introducing the replica symmetry breaking as 
pointed out in p. 

Another development of the theory should be the analysis of c? + 1- 
dimensional DPRM with d > 1. Some preliminary tests performed in d = 2 
pointed out that for the delta-correlated potential one must chose a proba- 
bility distribution which is not a simple power-law in k. Further inquiry are 
required maybe in the direction of a more complex functional structure of 
equilibrium probability distribution function. 

Acknowledgments I am deeply indebted with Professor Giorgio Parisi for 
his continuous support. 
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